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Paired state of nonstandard quasiparticles is analyzed in detail in two model situations. Namely, 
we consider the Cooper-pair bound state and the condensed phase of an almost localized Fermi 
liquid (ALFL) composed of quasiparticles in a narrow-band with the spin-dependent masses (SDM) 
and an effective field, both introduced earlier and induced by strong electronic correlations. Each of 
these novel characteristics are calculated in a self-consistent manner. We analyze the bound states 
as a function of Cooper-pair momentum |Q| in applied magnetic field in the strongly Pauli limiting 
case (i.e. when the orbital effects of applied magnetic field are disregarded). The spin-direction 
dependence of the effective mass makes the quasiparticles comprising Cooper pair spin distinguish- 
able in the quantum mechanical sense, whereas the condensed gas of pairs may still be regarded as 
composed of identical entities. The Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) condensed phase of 
moving pairs is by far more robust in the applied field for the case with spin-dependent masses than in 
the situation with equal masses of quasiparticles. Relative stability of the Bardeen-Cooper-Schrieffer 
(BCS) vs. FFLO phase is analyzed in detail on temperature - applied field plane. Although our 
calculations are carried out for a model situation, we can conclude that the spin-dependent masses 
should play an important role in stabilizing high-field low-temperature (HFLT) unconventional su- 
perconducting phases (FFLO being an instance) in systems such as CeCoIns, organic metals, and 
possibly others. 
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I. INTRODUCTION 



Unconventional superconductivity of heavy-fermion^ and organic metals^ is studied almost as intensively as that 
of high-temperature superconductors^.. While the rf-wave symmetry of the superconducting gap is shared among all 
of the systems, that have a quasi two-dimensional electronic structure, the first two systems are easier to approach 
theoretically, as the corresponding normal state can be conceptually described as a Fermi liquid, albeit an almost 
localized Fermi liquid, ALFL4, whereas the normal state of a high-temperature superconductor is an almost-localized 
non-Fermi quantum liquid, at least on the underdoped side^. By ALFL we understand the electron liquid located 
on phase diagram close to the Mott or Mott-Hubbard threshold, composed of (nonstandard) quasiparticles. The 
nonstandard characteristics of the quasiparticles comprise: (i) large value of their effective mass m* which becomes 
divergent if the Mott transition has the character of a quantum critical points, (ii) spin-direction dependence of 
the effective masses m* = m a in the magnetically polarized state in the non-half filled band situation 7 ., and (iii) the 
appearance of the effective field h cor induced by electronic correlations which differs from the molecular field introduced 
in magnetism 8 . Thus, it is of basic interest to include those novel features in the description of concrete physical 
phenomena, as they may mimic non-Fermi liquid features even though they represent only non-Landau corrections 
to the Landau Fermi liquid theory. The approach developed here is based on the Bardeen-Cooper-Schrieffer (BCS) 
theory of superconductivity and closely associated with it the problem of a single Cooper pair, both treated without 
and with the spin-dependent masses, in the latter case also with inclusion of the effective field h cor calculated within 
a self-consistent scheme. In this manner, our approach represents a natural and simple extension of both Landau 
Fermi liquid and BCS theories to the systems with correlated electrons. 

Spin- dependent masses (SDM) of quasiparticles have been observed recently in the heavy-fermion superconductor 
CeCoIn^ 9 . and other systems^ by means of the de Haas- van Alphen oscillations in strong applied magnetic field. These 
observations confirm the earlier theoretical prediction concerning the spin-dependent mass enhancement induced by 
inter-electron correlations in a single narrow-band 7 - and in the Anderson-lattice system s 11 ' 12 . This phenomenon has 
also been reinvestigated recently within the Periodic Anderson^ and the HubbarcU^ models. It appears that SDM 
should occur for either moderately or strongly-correlated systems with large on-site Coulomb repulsion U. This effect 
is particularly strong near the half filling n — ► 1 in the narrow-band or for an almost integer valency of Ce +4 ~ n/ 
(rif — > 1) in the heavy-fermion systems. Furthermore, in the same system - CeCoIn5, in which SDM were observed, 
a novel high-field low-temperature (HFLT) superconducting phase has been discovered^. This phase was proposed 
to be either of the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) natu re 16 ' 17 or an unconventional superconducting phase 
coexisting with an antiferromagnetic orde r 18 ' 19 , or even a complex phase with three independent order parameters^. 
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Hence, its nature is still unclear. Therefore, it is important to consider a case with SDM and reexamine the BCS phase 
against formation of the FFLO phase. Such analysis can help in predicting the influence of SDM on the HFLT stability. 
To achieve these goals we consider here only a model situation to emphasize the role of nonstandard quasiparticle 
characteristics in altering the standard BCS or FFLO description. In particular, we show that the FFLO state is 
robust in the SDM case on the expense of the BCS state. The simple model reflects some of the qualitative features 
of the phase diagram observed for CeCoIng^, although our main purpose here is to underline the universal properties 
rather than to construct the detailed phase diagram for concrete systems. In part of our numerical analysis, we take 
the values of some parameters (e.g. the elementary cell volume) such as for CeCoIn5 to illustrate that the results are 
in the proper range of temperature and applied field, although the detailed treatment for this system must be carried 
out separately and include also the d-wave symmetry, the singlet-triplet mixing in the paired state, and possibly, 
the spin-density-wave appearance, coexisting with the proposed here generalized FFLO state. We also discuss briefly 
some of the nontrivial features of the normal state. 

The structure of the paper is as follows. In Sec. II we analyze the concept of electron gas with SDM and of the 
effective field, and introduce its basic characteristics in the normal state. In Sections III and IV we discuss a single 
Cooper-pair problem^ for such Fermi liquid, particularly in applied magnetic field and obtain a stable pair state with 
center-of-mass momentum Q^O. The pair is discussed at some length, as it conveys a basic message concerning the 
formation of FFLO state and is nontrivial because quasiparticles composing it are spin distinguishable. In Sections V 
and VI we analyze the condensed state of pairs, both analytically and numerically. Finally, Sec. VII contains outlook 
and concluding remarks. In Appendices [A] - [C] we discuss the details concerning the origin of SDM, introduce the 
general antisymmetric state of Cooper pair, and justify the narrow-band limit for the case of Anderson lattice with 
an intraatomic hybridization, respectively. 



II. QUASIPARTICLE GAS WITH HEAVY SPIN-DEPENDENT MASSES AND EFFECTIVE FIELD 

INDUCED BY CORRELATIONS 

We analyze first the normal state properties of a quasiparticle gas with the spin-direction (rr = ±1) dependent 
masses m* = m a and the effective field induced by correlations h cor . Quasiparticle energies in the applied field 
h = g/isHa have the form 

% 2 k 2 

Ckcr = - — ah- (i- ah cor , (I) 

2m a 

where we have taken the simple parabolic dispersion relation and have defined from the start the energy with respect 
to the chemical potential \x. The spin dependence of the masse o 7 ' 12 is taken in the simplest form corresponding to the 
narrow-band or the Kondo-lattice limits (cf. Appendices lAl and [Cj) with the Hubbard interaction U — ► oo, i.e. 

m a l-n„ l-n/2 m _ 1 

= -: = -: cr;— - = (rriav -o-Aot/2), (2) 

tub 1 — n 1 — n 2(1 — n) tub 

where a = ±1 is the spin quantum number, is the band mass, rn = — n± is the system magnetic polarization 
and n is the total band filling (n = n-j- + nj,). Also, A771 = 777,2 — mi is the mass difference and m av = (mi + 7772 )/2 is 
the average mass. Note that in the magnetic saturation limit (77 j — n\)/ (n-j + rijj = 1 we recover the band limit with 
Tni/mB — 1, whereas the heavy quasiparticles in the spin- minority disappear (n.j_ = 0). Note also that the convention 
is such that the state a = +1 is regarded as that with magnetic moment along the applied field direction. 

In this paper we follow the spirit of original paper by Fulde and Ferrell^ and adopt it to the ALFL case. In 
such formulation, the system of self-consistent equations determining thermodynamic properties of the normal state 
starting from the free-energy functional J-, is as follows 

T = -k B T V ln(l + e- Kk °) +/iiV+ —mh cor , (3) 
z — ' 77 

kcr 
kcr 

™= Tf X^/^kcO, (5) 



n = 7i t = — ^/(£k<r), (6) 



N 
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FIG. 1: (Color online) Panel with the Fermi sea characteristics in the normal state as a function of applied magnetic field (for 
T — 0.05 if). Dashed line in a) represents the mass in the spin- minority subband, whereas the solid line characterizes that 
in the spin-majority subband. The dotted-dashed lines in d) represent the results for subband Fermi wave vector in the case 
with spin-independent masses (SIM) with m* = m av . Note much greater Fermi wavevector splitting in the SIM case; this is 
important for understanding of the results for superconducting state. For details see main text. 



where /(£k<r) is the Fermi-Dirac distribution, (3 = 1/ (fcgT) is the inverse temperature, and N is the total number of 
particles. The free energy functional T{/T, H a ; h cor ,m, n) given by (J3j) describes a Fermi sea with the spin-dependent 
masses m a and correlation field h cor . The equations (U) and (O are derived from the conditions dT /dm = and 
dT '/dh cor = respectively, and the last equation © is an elementary relation used to fix the band filling (defined by 
n/V e i em = N/V, where V e i em is the elementary cell volume). The normal-state properties determined via Eqs. ([3]) - 
(J6j) are to be compared with those for the paired state obtained in the subsequent sections. 

The equations describing the Fermi sea characteristics can be easily solved numerically by their reduction to a single 
equation for n„ of the following form 

njh^ njn-n,) 2 / 3 | 4(h + h cor ) / V \V3 

(n ~ n 1 )(2m av - m B ) + niBni ni(2m av - m B ) + m B (n - n\) h 2 \6ir 2 J 

with ni = jif. The Fermi sea characteristics are summarized in Fig. [TJ The mass difference, the Fermi vector 
splitting and magnetization increase linearly with the increasing field. Therefore, the approximated expressions in 
the field up to 30 T are m = xH a and m a {H a ) = m av — <jmBx/2(\ — n) H a . In our numerical study we obviously use 
full expressions ([3]) - ([6]). Namely, we take the values of parameters emulating the heavy-fermion systems: n = 0.97, 
V e i em = 161 A 3 and the h = value of the quasiparticle mass m av = raj 1 ~"^ 2 = 100 ttiq. For these values the principal 
characteristics are collected in Fig. Q] The assumed mass enhancement magnitude corresponds to the moderate heavy 
fermions, with the value of 7 in the range 100 — 200 mJ/mol K 2 . Also, the value of n = 0.97 corresponds to the 
effective valence of the Ce ions +4 — nf = 3.03, a typical value. Note that the mass splitting is only about 7% in the 
field of H a = 30 T, but more important is the Fermi wavevector splitting Afcp = kp-\ — kpi, displayed in Fig. [TJi. 
The Fermi wavevector is calculated according to the relation kp a = (67r 2 n cr /V r e / em ) 1 / 3 . Most of the characteristics are 
indeed linear in H ai as stated above. 
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FIG. 2: Correlation field h cor as a function of the applied magnetic field for the normal state. The linear dependence is 

hcor = —0.4 A. 

In Fig. [5]we display the h dependence of the effective field h cor . It is linear in h and typical values are h cor w —0.4 h. 
More importantly, it is always antiparallel to the applied field, it partly compensates it in the sense that it reduces 
the Zeeman contribution to the quasiparticle energy. Also, the external field induces the effective-mass splitting and 
this factor drastically decreases the difference 

III. COOPER PAIR STATE: THE QUASIPARTICLE DISTINGUISHABILITY 

We consider next the simplest situation involving quasiparticle pairing - the Cooper problem 2 ^ and demonstrate 
principal and nontrivial features of the introduced quasiparticle states, particularly in the applied magnetic field. 
These features involve not only the pair binding energy, but raise also some basic questions concerning its quantum 
mechanical description, as discussed next. 

A. Hamiltonian with scalar effective masses and its limitations 

A simple way of accounting for the different masses of quasiparticles is to insert them into the Hamiltonian "by 
hand" , that is to assume that the first particle has the spin quantum number a — \ = 1, and the corresponding mass 
m-f = mi, whereas the other one has spin a =1 = — 1, and mass mi = 7712. This is a standard procedure, since in 
nonrelativistic quantum mechanics the particle mass is an external, classical parameter. This assumption leads to the 
following form of the Hamiltonian for the two particles with m-\ = mi and = mi 



H = V? - ^— V 2 2 - afh - of h - afh cor - a z 2 h cor + V(n,r 2 ). (8) 

Where V r (ri,r2) is an attractive interaction between them. We neglect the effect of applied magnetic field on the 
electron orbital, because the Maki parameter— is usually large in correlated systems (see e.g. CeCoIns, in which 
a" = 4.6, a 1 - = 5.0i^), which means that the Pauli contribution dominates over the orbital effects. For the two- 
particle state with opposite spins, Hamiltonian ||5J) acts properly only on the wave function with the spin part of the 
form 

x(oi,o 2 ) = |lT)|2i>, (9) 

This type of wave function will be called specific-spin in the following. If we applied the Hamiltonian to the wave 
function with different spin part, i.e. x'( <T ij (J 2) = |1 l)|2 T); it would assign improperly mass mi = m-\ to the first 
quasiparticle, whose spin is a =j. For this reason, we may not construct for Hamiltonian ([5]) the proper singlet or 
triplet spin wave functions Xs,t(oi, 02) = (|1 t)|2 |) T |1 l)|2 T))/\/2, in which the spin quantum number has not 
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been assigned to the particle of given mass. This means that by specifying that the first particle has the spin up, and 
the second the spin down we violate in an obvious manner the quasiparticle-spin indistinguishability. 

One can illustrate the nontrivial character of the present case with SDM by commenting on the simplest quantum 
system - the hydrogen atom. Namely, although the masses of proton and electron are vastly different, the total spin 
part of the wave function is still the full singlet. This is because the spin quantum numbers of those particles are 
completely detached from their masses. In effect, the spin part of the wave function is either pure singlet or triplet. On 
the contrary, in our situation with SDM, the spin quantum number is attached to the masses, so the spin transposition 
symmetry may be broken explicitly, as discussed in detail below. 



B. Hamiltonian with masses in an invariant (operator) form 



We construct next the Hamiltonian in such a way, that it properly assigns masses to quasiparticles depending 
on their spin direction in the applied magnetic field. The only way to do this is to introduce the mass operator 
rh(af) = m av — ^o~f Am, where the mass splitting Am needs not to be specified at this point. Under this prescription, 
the two-particle Hamiltonian takes the form 



H = ~n z~\ ^1 - o TT—^l - °l h - a 2 h - aZ l h cor ~ O^hecr + V(n,r 2 ). (10) 

2m av — erf Am 2m av — erf Am 

Now, the kinetic part assigns respective masses mj = m av — Am/2 and m 2 = m a „ + Am/2 to the particles depending 
on their spin z-component. With this Hamiltonian we can now analyze any spin function. In particular, both the 
specific-spin wave function ([9]), as well as those describing the singlet or triplet states with the z-component of the 
total spin Sf + Sf = 0, which have obviously the form 



XS,t(°U°3) = ^(l 1 t)|2 I) T |1 i)|2 T>). (11) 

Additional advantage of the Hamiltonian (JTDJ) is the ability to describe transition from the indistinguishable- 
quasiparticles limit (for zero or low mass splitting) to the spin- distinguishable- quasiparticles case when h > 0. By 
distinguishable quasiparticles we mean here those described by a wave function without well-defined transposition 
symmetry (i.e. non-antisymmetric in our case). 



C. Singlet, triplet, and their inadequacy in describing the paired state 

First, we note that the total spin operator S 2 = (Si + S2) 2 does not commute with Hamiltonian (JTDJ), i.e. 



[H, (Si + S 2 ) 2 ] 0, (12) 

whereas H. does commute with Sf and separately. This means that while the z-components of the individual 
spins represent good quantum numbers, the total length does not. To analyze this property, it is important to see 
that the spin-dependent denominators can be rewritten in the two equivalent forms 



n = J2 

i=l 
2 

= E 



o i~r~~~V 2 - v*h - o--h cor 

2m„„ — afAm 



r 2 



+ V(ri,r a ) 



4m 2 )J — (Am)' 



-(2m av + erf Am) V 2 -afh- <r*h 



+ V(ri,r 2 ). 



(13) 



Taking Si = (l/2)<fi, one can easily prove the condition (TT2"]) . as well as the property that [H, Sf] = 0. One sees that 
the spin wave function can be characterized by individual values a\ =\ and erf =1 or vice versa, but the two-particle 
spin state might not have a proper singlet or triplet symmetry (111|) . Note also that the property (|12p is independent of 
the form of pairing potential and is fulfilled also for V(r 1 , r 2 ) = 0. This may lead also to the normal-state corrections 
not discussed here. 
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D. Solution for the specific-spin wave function 

The wave function is decomposed into the spin and space parts 

*(ri,r2,CTi,cr 2 ) = $(ri,r 2 )x(o-i,c r 2 ), (14) 
with the space part given as a superposition of plane wave states, i.e. 

$(ri,r 2 ) = £ ^^^(rO^Cr,) = ± £ a kl , ka e^n+to (15) 

ki,k 2 ki,k 3 

By starting from the wave function in the form (114[) with the space part given by (|15|l and spin part x(<Ji, cr 2 ) = 
|1 |)| 2 |), we can solve the Schrodinger equation (with any of the two Hamiltonians) in a similar manner as in the 
original Cooper problem^.. One of the differences we introduce is a new definition of the relative momentum when 
transforming the Hamiltonian to the center-of-mass and relative coordinates, namely 

rjmi + r 2 m 2 

R= ■ ; Q = ki+k 2 , (16) 

7711 + 7772 

kl777 2 -k 2 mi 

r = ri-r 2 ; k= ■ , (17 

777l + «2 

where Q is the total momentum and k is called, by analogy to the standard Cooper problem, the relative 
momentum 23 . After this transformation the Hamiltonian and wave function can be cast to the forms 

H = -^rrV^ - ^— V 2 - a{h - a z 2 h - a z x h cor - a z 2 h cor + V{r), (18) 

2M Zfl re d 

r, oi, er 2 ) - -^e 4QR ^a k e lkr |l T>|2 I), (19) 

k 

where M = mi + 777 2 and /J, re d — 777im 2 /M. Following the standard procedure, we obtain the equation determining 
coefficients a k and the eigenenergy E in the form 

<»> 

where £q = ?i 2 Q 2 /2M , e k = ?7 2 k 2 /2/j re( i and N is the total number of particles. The interaction region has to 
be defined in a more general way, since we want to describe a system with nonequal Fermi vectors and possibly, 
with a non-zero center-of-mass momentum Q, which is a constant of motion. We assume that a constant, attractive 
interaction takes place in the regions of k-space for which both particles are at most at the distance htoc above their 
Fermi surface, i.e. 



[0, in other cases, 



(21) 



where k and Q are related to ki, and k 2 via transformation (fT6|) and (fl~7|) (the same holds for the vectors pairs 
k', Q and k' l5 k' 2 ). We call Wq the region in k-space in which the interaction is nonzero (cf. Fig. [3]). It can be shown 
that Wq — [Wi n (W 2 + Q)] — ^m"Qj where by adding a vector to the region in k-space we mean the whole region 
shifted by that vector. Also, Wi — {k | < £ k(Ti < hue}, i = 1,2. In this notation, equation for the binding energy 
A = 2/i — E becomes 

N V f d 3 k V f d 3 k 

(22) 



V 8^ J Wq e Q + e k - 2 M + A 8tH» J Wq Ck lT + Ck ai + A 



FIG. 3: (Color online) Interaction region Wq is the common part of W\ and W% + Q shifted by — ^Q. The vectors k a i, are 
defined by £k al t — £k a2 -L — ^c- Regions which contribute most to the pairing are marked as the three (green) ovals on the 
right with £,k iCri « 0. For details see main text. 



From this form of the equation for A we can deduce that the regions of reciprocal space contributing most to the 
pairing are those for which £k i(Ti — (see Fig. [3]). Large part of the space fulfills this condition if |Q| ~ Akp = 
kpi — kp 2 - Therefore, we can anticipate that the pair will have maximum binding energy when the pair center-of-mass 
momentum is close to the Fermi vector splitting, i.e. when |Q| — Akp, consistent with the results for the FFLO state. 

The equation (|22j) for the gap A has to be solved numerically for each |Q| and the final solution in our approximation 
is the one with the largest binding energy. In the case of Cooper pair at rest (Q = 0), equation (|2"2")l can be solved 
analytically. Such analytic solution depends on which of the two vectors k a i = ^2mj(/i + Oih + <Jih CO r + hcuc) is 
larger, namely 



A = . K h + h *or)M ~ [lAm] -^- l ( < ±- + J-e c ) } for k al > k a2l ^ 

ft "e_T - ^ [(h + h cor )M - fiAm] , for k al < k a2 , 

where C = v - ■ In zero magnetic field we have that C — yrjj—, and expression (|2"5]) properly reproduces 

the Cooper result A = 2hu>c / [exp( y p^ eF ) ) — l]j where ~p(ep) is the average density of states per particle. 



IV. DETAILED ANALYSIS OF THE COOPER-PAIR STATE 



A. Is the singlet state a proper eigenstate? 

As has already been said, Hamiltonian (fTU|) does not commute with the total spin S 2 . This means that the singlet 
and the triplet states are not a good basis for the problem considered. Nevertheless, we can still analyze energy for 
the singlet spin function (and we can think of it as a superposition of eigenfunctions with spin parts |1 |)|2 |) and 
|1 j,)|2 I)). This means, we can take a particular solution in the form 



*(r 1 ,r 2 ,a 1 ,a 2 ) = l £ « kl)k2 e *^+*a*» -L (|1 T )|2 |) - |1 |)|2 |». (24) 

ki,k 2 

Solution of the Schrodinger equation in this case is quite cumbersome and will not be presented in detail here. It 
can be shown that, for the resulting equations not to be contradicting, the interaction has to be introduced as nonzero 
in the region (for derivation see Appendix |B|) 



Cooper pair at rest 

exchange 
of spins 



4 i i 4 

00© - ©(2) 

moving Cooper pair 
Q#0 



exchange 
of spins 



CO© " ©(2) 



FIG. 4: (Color online) Behavior of the spin part of the wave function under transposition of particles labeled 1 and 2 for Cooper 
pair at rest (upper row, antisymmetry) and moving Cooper pair (lower row, no definite symmetry). 



w * = ^ = kima — kami I ^ e n W2 A k 2 € W2 H Wi A Q = ki + k 2 }, 

mi + ni2 



(25) 



or equivalently, 



W* 



w Q n(-w Q + Q- 



(26) 



The replacement of Wq with W* presents the only difference in the case of pure singlet state, with respect to 
the that for the specific-spin state |1 f)|2 J.) or |1 J,)|2 |). This interaction region coincides with Wq if and only if 
Q = 0. As a result, the energy of the singlet state in the case of Cooper pair at rest is identical to the energy for 
the specific-spin states. For a moving pair (Q 7^ 0) the interaction region is reduced, and the binding energy for the 
singlet state decreases rapidly with increasing |Q|, whereas that for the specific- spin states can become even larger 
(see Fig. |3J}. 

B. Quasiparticle distinguishability 

The wave function of fcrmions has to be antisymmetric with respect to transposition of spin and space coordinates 



(ri,o-iJ <-> (r 2 ,cr 2 ) 



(27) 



For the Cooper pair this implies, that 



*(ri,r 2 , 01,0-2) = -*(r 2 ,ri,cr 2 ,o-i). 



(28) 



It turns out that the Cooper pair at rest (with the singlet spin part) has proper symmetry even for mi 7^ rri2, a 
very interesting result. As long as Q = leads to the minimum energy state, the wave function is antisymmetric and 
describes indistinguishable quasiparticles. 

For sufficiently strong applied magnetic field the state with Q / becomes stable and the pair has to be either 
in the |1 |)|2 j.) or |1 J,)|2 f) state. Such wave function does not have the symmetry (f2"5]) and the transition to 
spin- distinguishable quasiparticles takes place. The difference between the two situations has been illustrated in Fig. 
S] and the former providing the connection of the wave- function symmetry with the type of the solution. 

A word of explanation is most proper at this point. Namely, the microscopic many-body Hamiltonian (Hubbard 
or Anderson-lattice, for example) respects the particle indistinguishability. This microscopic model is approximated 
here by an effective quasiparticle picture with the SDM£ i 12 ' 13 i 14 . The effective Hamiltonian taking into account 
those masses violates the indistinguishability in an obvious manner (see relation (fl"2"]) ). However, this quasiparticle 
picture has been confirmed experimentalljs^i for the normal state. Therefore, this fact suggests that there is a basic 
qualitative difference between the Landau Fermi liquid describing rather weakly interacting fermions and the almost 
localized (local) Fermi liquid more appropriate for moderately or strongly correlated particles. 
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FIG. 5: (Color online) Binding energy of the specific-spin (solid lines) and singlet states (dot-dashed lines) as a function of 
center-of-mass momentum Q and for selected values of applied field. Note a rapid decrease of binding energy of the singlet 
state. The bound state with |Q| 7^ becomes stable above the field H a — 2 T. The arrows mark the state with maximum value 
of binding energy. 



C. Numerical results 



In numerical analysis of the Cooper pair state we assume the following values of parameters, emulating the heavy- 
fermion systems: n — 0.97, V e iem = 161 A 3 , m av = ms \^ = 100 m (all taken for CeCoIn 5 ). Also, we take 
Tiljc = 100 K and Vb = 93 X. Additionally, for sizeably different values of hue, another maximum in the binding 
energy appears (between Q = and |Q| = Akp), with binding energy maximally 5 mK higher than that for Q = 0. 
Therefore, with this third solution no new physics is incorporated in our study, and it might only blur the image 
presented. For the assumed parameters, the chemical potential is /i rs IA0K, i.e. of the same magnitude as hue and 
Vq. The relative value of the parameters is of no primary importance here, since the solution is of nonperturbational 
nature. 

Solution of the gap equation (|22|) provides us with dependence of A on the Cooper pair center-of-mass momentum 
Q, in an applied field H a , as shown in Fig. [5] For the fields above 2T the maximum binding energy appears for 
|Q| ~ Akp. Note that the full singlet configuration (dot-dashed lines) has much smaller binding energy than the 
specific-spin state (solid lines). 

The binding energy as a function of the field is plotted in Fig. |6] The case with SDM gives rise to much higher 
critical fields above which the pair is destabilized by the Pauli effect. The reason behind this robustness of the state 
for Am ^ is the smaller Fermi wavevector splitting Akp for this case (c.f. Fig. [l}, which in turn leads to the larger 
interaction region Wq. One can say, that the effect of SDM acts in the opposite direction than the Zeeman-term 
influence and compensates this influence to a degree. 

The optimal momentum value |Q| vs. H a is displayed in Fig. [7] One observes a discontinuous change of |Q| 
above the field ~ 2T from |Q| = to |Q| ~ Akp. Above this point the specific-spin wave function provides the 
stable solution. At the same time, we have a transformation from the spin-indistinguishable to spin- distinguishable 
quasiparticles. In this manner, we have a model situation in which the question of distinguishability of quasiparticles 
can be investigated. 



D. An overview of the Cooper pair problem 



We have presented solution of the Cooper pair problem with SDM. The Hamiltonian has been formulated in two 
ways: with the scalar masses and with operator masses, the latter being more appropriate for the study of the 
two-particle pairing. The latter Hamiltonian also allows for systematic studies of the question of quasiparticle-spin 
distinguishability. The resulting stable solution has an antisymmetric wave function and describes spin-singlet state 
for Cooper pair at rest (Q = 0), even if mi ^ m-2. However, for (and only for) Q = the singlet solution is 
energetically degenerate with the specific-spin solution. For higher magnetic field, the solution with Q = becomes 
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FIG. 6: (Color online) Pair binding energy as a function of magnetic field for the cases with the spin-dependent (solid lines) 
and the spin-independent (dot-dashed lines) masses of quasiparticles. Note higher critical fields (as marked by the solid circles 
on the s-axis) for pair breaking in the former case. Analytic solution reproduces properly the dependence for Q = state. The 
uppermost solid circle marks the transition from Q = to |Q| ~ Afc^. The state with Am — and Q 7^ is nonphysical and 
is presented for comparison. 
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FIG. 7: (Color online) Optimal Cooper pair center-of-mass momentum versus magnetic field. For high enough fields |Q| ~ AftF- 
The regimes with spin-singlet and specific-spin wave functions are also shown. The solid circle at the end marks the Pauli 
limiting critical field. 



energetically unfavorable, and this is why the transition to a moving-Cooper pair takes place. The moving Cooper-pair 
wave function has no definite transposition symmetry in the spin coordinates. Therefore, this transition is a typical 
transition to a broken symmetry state in the Landau sense, since the starting Hamiltonian with the effective masses in 
the operator form (cf. Eq. (|10[) ) is written in spin-symmetric form. With the transition to a moving Cooper-pair state 
there also appears that from indistinguishable quasiparticles to their spin- distinguishable correspondants. This also 
leads to the admixture of the triplet component with S z = to the singlet state for Q ^ 0, as the specific-spin function 
is a superposition: |1 f)|2 J.) = Qsinglet) + \triplei)) / V2- We think that the quasiparticle-spin distinguishability is 
thus indispensable, as the quasiparticles have different masses which represent their external characteristic, with the 
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value of spin quantum number labeling them. 



V. PAIRED STATE OF QUASIPARTICLES IN ALMOST LOCALIZED FERMI LIQUID 



A. General remarks on the pairing nature 



At present, two principal non-phonon mechanisms of pairing have been applied to strongly correlated electron 
systems. The first of them is the real-space pairing based on either t-J model and applied to high temperature 
superconductivity 2 ^ or on hybrid pairing applied to either heavy- fermion— or high-Tc superconductivity 26 . The 
second is based on (renormalized) paramagnon exchange 27 . In both cases, the effective pairing potential is explicitly 
momentum-dependent and this situation results in k-dependent gap. For the purpose of the present paper, in which 
we take the narrow-band limit of the Anderson lattice model with the k-independent hybridization matrix element (cf. 
Appendix [C| , the pairing potential near the Fermi level is approximately constant and thus leads to the s-wave form 
of the superconducting gap. Note however, that the original purpose of the present work is to formulate the language 
of description of the systems with SDM and h cor present in a simple situation. Such straightforward formulation leads 
already to quite nontrivial and probably universal mechanism of the HFLT phase stabilization by SDM (FFLO being 
an instance). 



B. Formulation of the problem 



We discuss next the condensed state of the nonstandard quasiparticles introduced in Sec. [TTJ In the paired state 
with isotropic gap the applied field penetrates only to the London depth A^, i.e. produces an irrelevant contribution 
of the spin-split-mass states onto bulk properties at temperature T = 0. However, at T > the excited quasiparticles 
will produce a nonzero moment which should be clearly visible when the system approaches the transition temperature 
T a . Moreover, SDM should provide an important contribution to the FFLO superconducting state, since in that phase 
there are normal regions. Therefore, in the situation with SDM (and with h cor reducing h) we should observe a robust 
FFLO phase in high fields on expense of the BCS state, as SDM help stabilizing it (by reducing the Fermi vectors 
splitting AltF and by increasing the critical field due to the presence of h cor ). This is the subject of discussion in the 
next two Sections, as a result of which we construct the phase diagram on the applied field H a - temperature T plane. 

We introduce the BCS Hamiltonian with a constant pairing potential of magnitude Vq (for its simple justification 
see Appendix [C]) and allow for the possibility of a nonzero center-of-mass momentum Q of a Cooper pair written in 
the representation symmetric with respect to both component quasiparticles composing the Cooper pair— 



W — XI € k <r a la ak * ~ 77 X a k+Q/2T a -k+Q/2l a -k' + Q/2|ak' + Q/2T + —mh CO r, (29) 
kcr kk'Q 

with £ka- given by (fl]) taking account of the different quasiparticle masses. The magnetic field is accounted for 
only via the Zeeman term, as the Maki parameter— in the systems of interest is highi^ (Pauli limiting case). The 
interaction is assumed to exist only in region ±hu>c around the Fermi surface, more precisely in the region 



_ = - fc bT + k bl k a] + k al - ^ 

where kb a is defined by £k ba ,cr = —fvuci an d k a( j, by ^k atT ,a = hu>c- Such interaction region has been chosen because 
its width does not change significantly with the magnetic field h. We performed also calculations by selecting the 
interaction regime differently (namely, by choosing W = [fc&^, fe a f] and W = [k^, fc a jj) and have obtained almost the 
same results. Hamiltonian (f!?9")) is diagonalized with standard mean-field procedure followed by the Bogolyubov-de 
Gennes transformation of the form 

Q^kT = "kak+Q/2T ~ w kal k+ q/2X; (31) 
at t = w k a k+ Q/2T + M ka^ k+Q/2i , (32) 



which leads to the diagonal form 
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H = E^al^a + ^(Ck } - Ek) + N-r^- + -mh cor , (33) 



and the quasiparticle spectrum characterized by energie: 



,28 



^k 
1 



e k s)2 



A 2 



= ^(£k+Q/2T + £-k+Q/2J.), £ k ° = ^(£k+Q/2"T ~ £-k+Q/2J.)> 



(34) 
(35) 



and with 



w2> 



■k+Q/2|ak+Q/2T) 



(36) 



being the Q-dependent gap parameter. In the form l|31 [l - (|32p of the Bogolyubov-de Gennes transformation the 
quasiparticle operators a k f and a^, are distinguished by the pseudospin label f and J,. Note also that because of 

the presence of £ k , there are regions of reciprocal space for which _E k(T < 0, which represents nongapped excitations. 
The gap parameter Aq is determined from the self-consistent gap equation 



1 - /(£ kT ) - /(gkj) 

2E k 



(37) 



The Bogolyubov transformation coherence factors in (l3Tj) and (|32|) are given by 



— ) 



1/2 



f'k 



—) 



1 1/2 



(38) 



Finally, the complete set of equations determining the superconducting state properties is as follows, 



T = -fc B T^ln(l + e -^)+^(^. 

kcr k 



•JV- 



Vn 



iV 



kcr 



dm N ' 9to 

k 



m 



kcr 



vb^ i-/(£; kT )-/(g ki ) 
^ 



2£ k 



TV 



^{ Wk /(£; k(T ) + ^[i-/(^ k) - CT )]} ; 

kcr 



(39) 

(40) 
(41) 

(42) 
(43) 



where T(T, H a ; h cor ,m, n, Aq) is the system free energy for the case of a fixed number of particles^. Similarly as 
for the non-interacting Fermi sea, the equations ([40]) . (|4Tj) and (f42|) are equivalent with dT '/dm — 0, dJ-/dh cor = 
and 8J-/dAQ = 0, respectively. In effect, the numerical analysis involves solving the system of four integral equations. 
Note also the presence of two different effective chemical potentials fi a = fi + uh cor for particles with spin up and 
down in the spin-polarized situation. This is an unavoidable consequence (c.f. Refsf£ & 12 i 30 ) of the slave-boson 
formalism used to derive expression for the masses @, and dispersion relation |T]). Parenthetically, we have also 
performed calculations by disregarding the different effective chemical potentials (i.e. we have put h cor — 0) and the 
results obtained were nonphysical (the free-energy jump occurred at the BCS-FFLO phase transition). Obviously, 
the physical chemical potential is [i and is determined from the neutrality condition (|43p . 

The physical solution is that with a particular Q which minimizes the free energy (|39p . The state with Q = is 
called the BCS state, and that with |Q| ^ - the FFLO state. 
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Phase Diagram - SIM 
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FIG. 8: (Color online) Phase diagram for the cases with the spin-dependent a) and the spin- independent masses b). Light 
(yellow) region corresponds to Q = (BCS phase), the darker (blue-red) one to Q / (FFLO phase) and the white to normal 
state. Note that with increasing temperature, the transition from BCS to FFLO state occurs at higher fields, in qualitative 
agreement with experimental results-^. The FFLO phase is stable in an extended H a -T regime only in the SDM case. 



VI. NUMERICAL ANALYSIS AND DISCUSSION 



We assumed, the following values of the parameters, emulating the heavy fermion systems: n = 0.97, V = 161 A 3 , 
m av — 100 mo (data for CeCoIns), %ujc = 17 if, and Vo = 93 if. The characteristic energy scale associated with 
spin-fluctuations in CeCoIns is T s f = 10 if— - a value comparable to our hu>c- For those parameters, the chemical 
potential was equal to /i « 140 if. This means that Vo ~ ep and the (weak coupling) BCS approximation can be 
regarded only as a proper solution on a semiquantitive level at best. Furthermore, for the values of parameters one 
can calculate the coupling constant p(fi)Vo ~ 0.48 and the coherence length at T = if , £o ~ 40A, both already at the 
border of the strong-coupling limit. Also, such values have been taken to obtain the critical temperature T s ~ 2.3if . 

We now discuss the phase diagrams for the cases of SDM and SIM which are exhibited in Figs. [SJi-b, respectively. 
Both the BCS (state with Q = 0) and the FFLO (Q ^ 0) phases extend to much higher fields if the masses are 
spin-dependent. This is a consequence of the smaller Fermi-vector splitting Akp for the SDM case (c.f. Fig. [TJi) . 
The most interesting is the fact that in the SDM situation the FFLO state becomes much more robust compared to 
BCS state, especially for T ~ 0. The reason for this is as follows: the superconductivity in the Pauli limiting case is 
destroyed by the Fermi vectors splitting Akp (c.f. Fig. [TJi). This splitting in the case of SDM is generally smaller (in 
this respect SDM compensate effect of the Zeeman term), hence the higher critical fields. However, for the masses 
to depend on spin, the magnetization m = nj — nj, has to be non-zero, and in the BCS state around T < 0.5 if 
magnetization is close to zero (see Fig. [5]) what weakens the mass dependence on spin and in effect produces larger 
AkF- Therefore, the BCS state is not enhanced much by SDM in that temperature interval. In the FFLO state, on 
the other hand the magnetization is nonzero even at T = Oif. This is because in the FFLO state there are regions 
with unpaired quasiparticles in the reciprocal space. The FFLO state becomes stable in this regime, as a result of a 
smaller Fermi vectors splitting. 

Another interesting feature is the fact that with increasing temperature, the transition from BCS to FFLO state 
occurs at high fields (cf. Fig. EK) consistent with experimental results^ 5 . It can also be easily explained. As 
temperature increases, the magnetization in the BCS state increases, allowing a substantial mass difference, and 
decreasing the Fermi wavevectors splitting, enhancing superconductivity. Therefore, the BCS state benefits from the 
smaller Afejr for SDM at higher temperatures (T > 0.5 if) and becomes more stable in this regime. 
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FIG. 9: (Color online) Spin magnetization as a function of temperature and applied magnetic field. For T = — 0.5 K, the 
magnetization in the BCS state is small, than it increases with applied field after the continuous BCS to FFLO transition. 
For higher temperatures magnetization in the BCS state becomes substantial and this produces a higher critical field for the 
BCS-FFLO transition for T > 0.5 K. 



Systematic evolution of the pure spin magnetization in the condensed state is shown in Fig. [9] (the orbital part is 
not included). It increases at the BCS-FFLO border al lower T, as one would expect. 

In the panel composing Fig. [10] we plot the gap magnitude Aq and the magnitude of the wave vector Q for SDM 
and SIM cases. The behavior of the order parameter Aq differs substantially in these two cases. Namely, there is no 
jump of Aq at BCS-FFLO transition for SDM, whereas for SIM this transition is always discontinuous. Transitions 
from superconducting to normal state are continuous for the case of SDM in disagreement with the experimental 
results^. The reasons for this discrepancy are discussed in the next Section. 

Finally in Fig. [TT]we show the correlation- field dependence. It can be seen that for BCS around T = OK this field 
is close to zero, then increases and approaches for H a — > H C 2 the value for the unpaired Fermi sea, denoted here as 

hcorFS- 



VII. OUTLOOK AND CONCLUSIONS 



We have analyzed superconducting states of a three-dimensional gas of heavy quasiparticles with (SDM) and 
without (SIM) the spin-dependent masses. Despite the simplicity of our model (parabolic dispersion relation, constant 
pairing potential, s-wave gap, single narrow-band model), qualitative results obtained are very meaningful for the 
FFLO phase detectability and should hold for more general and realistic models and other high-field low-temperature 
(HFLT) phases. This is because the spin-dependent factor renormalizing mass is k-independent, as is h cor , and they 
are obtained in a self-consistent manner from global equations, integrated over k. In effect, their values should not 
be influenced strongly by the details of the bare electronic structure. 

The most striking result is the fact that for the case of SDM the FFLO state becomes stable in much wider 
range of applied field and temperature. We believe that the mechanism of stabilization of the FFLO state by SDM 
is universal. Therefore, it should also apply to other unconventional HFLT phases such as for example the mixed 
staggered 7r-triplet SC + d-wave singlet SC + SDW phase proposed very recently 2 ^. This is because in such phases 
the spin-magnetization is always higher than for the conventional BCS (s-wave or d-wave) stat o 20 i 32 . Those phases 
will benefit, even to a larger extent, from the compensation of Zeeman effect by SDM, on the expense of the BCS 
phase, as discussed earlier. 

The detailed application of our results to concrete systems is rather limited. This is because three topics require 
still a conjoined analysis from the theoretical side. This is the inclusion of the singlet-triplet mixing in the FFLO 
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FIG. 10: (Color online) Left (a, c): gap parameter Aq as a function of temperature and magnetic field for SDM (up) and 
SIM (down). All transitions for SDM are continuous. Right (b, d): Cooper pair momentum in the FFLO state in units of 
Fermi- wavevector splitting AIcf for SDM (up) and SIM (down). Note that for the FFLO phase the momentum |Q| changes 
continuously (with transition BCS-FFLO), contrary to the case of SIM. Typical value of the momentum is |Q| ~ A/cf- 

phase?— as mentioned in the Cooper-pair case (cf. Sec. IIIip . Associated with it is the problem of magnetism 
appearing in the normal portion of the system in the FFLO state. Third, we have to introduce ti-wave symmetry of 
the superconducting gap. Inclusion of those factors introduces additional self-consistent integral equations making 
the whole approach much more complex from the numerical side. Then, one has also to carry out the whole procedure 
for a realistic (quasi-two-dimensional) electronic structure. The inclusion of magnetism should result in the first-order 
nature of the BCS-FFLO phase boundar y 19 ? 20 . We should be able to see a progress along these lines in the near 
future. 

The nature of the HFLT unconventional phase in the heavy - fermion system CeCoIns is still unclear. Some studies 
suggest FFLO character—, others reject it^. We claim that whatever this state really is, it may be stabilized by SDM 
due to its higher spin-susceptibility 32 . So far the observation of FFLO phase in organic metal has been confirmed 3 ^, 
but no spin-dependence of the effective mass has been investigated for those systems. 

To conclude, the simultaneous observation of the spin-dependent masses and of an unconventional HFLT supercon- 
ducting phase in the same system should not be regarded as coincidental. Hence, other unconventional HFLT phases 
may be searched for in the systems in which spin-split masses have been observed and vice versa. 
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FIG. 11: (Color online) Left: Correlation field in units of the field obtained for non-interacting Fermi sea with SDM included. 
Note that as H a — * H C 2, the correlation- field value approaches the one for the unpaired Fermi sea. Right: Correlation field in 
absolute units. It is negative, i.e. it acts opposite to the applied magnetic field. 



APPENDIX A: ELEMENTARY DERIVATION OF THE SPIN-DEPENDENT MASS ENHANCEMENT 

FACTOR 

The mass enhancement factor in a strongly correlated system can be understood on the intuitive ground in the 



.7.13.14 



, so 



following manner. Usually, it is the most decisive factor in determining the quasiparticle density of state: 
the argument is carried out for featureless (rectangular) form of the density of states in the bare band (here taken 
per site and one spin direction), i.e. 



f J_ for - K < f < K 
p(e) = <^ w 1 2 - - 2 

ry ' [ 0, otherwise, 

where W is the bare bandwidth. In such situation the bare band energy per site is 



(Al) 



V N Ju^q 

where the chemical potential is defined from 



ep(e)de 



(y) *52 n M ~ n*)> 



(A2) 



p(e)de. (A3) 
On the other hand, for strongly correlated electrons (U — > oo) the corresponding energy of itinerant electrons is 



/Eb\ (W 



\ N Ju^c 



(f)E»^-»)' 



(A4) 
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The factor (1 — n) expresses the fact that the hopping of the electron between the neighbors takes place when there 
is no other electron present on neighboring site. Combining (|A2p and (|A4p one can write down that 



(^)^oo =_ (T)^^ (1_?lff)= ^ ^ £k ' (A5) 



h<ihp ,er 

with q a = (1 — n)/(l — n a ) is the spin dependent, k-independent renormalization factor. We write now the 
quasiparticle energy in the form = q CT £k = + {q a — l)ek = £k + S CT (ek), where S CT (ek) is the self-energy part 
(its real part) induced by the correlations. Defining the mass via the standard Fermi-liquid relation 

+ (A6) 

we obtain that m* = m a = ms and \/q a — to_b(1 ~ ?v)/(l ~ which is equivalent to The spin-dependent 
masses of this form appear in either Gutzwiller— or slave-boson approaches^ 2 - which reduce to the Gutzwiller-ansatz 
results in the paramagnetic case. They do not appear in the simple slave-boson approach^, but then the overall 
quasiparticle mass depends on applied field. The appearance of the effective field h C07 £ & 12 i 30 arises from a constraint 
on the number n a of fermions, in direct analogy to the introduction of chemical potential in grand-canonical approach. 

APPENDIX B: PAIR BINDING ENERGY FOR GENERAL ANTISYMMETRIC STATE 

We consider here the most general wave function antisymmetric with respect to transposition of particles (|27[) and 
show that its binding energy decreases with the increasing center-of-mass momentum Q and is lower than the energy 
of the specific-spin state ©. In the main text we refer to this solution as to the singlet (which is a specific case of 
the considered here general function), because we would like to recover the standard singlet solution in the H a = 
limit. Spin states with a well-defined transposition symmetry (singlet and triplet with S z = 0) are described by the 
following wave functions 

Xs(a 1; a 2 ) = i=(|lt)|2|)-|U)|2T)), (Bl) 
X T(a 1 ,a 2 ) = -i=(|lT)|2|) + |U)|2t)). (B2) 
Therefore, we can expand the general antisymmetric wave function in the basis of singlet and triplet wave functions. 



\£r(ri,r 2 , (Ti,(T2) = A$5(ri,r 2 )x5(c r i,c r 2 ) + yl - A 2 <&a{yi, r 2 )xr(<7i, a 2 ). (B3) 

Where $5 is symmetric and $,4 antisymmetric under (ri <-> r 2 ) and A is characterizing the degree of mixing. We 
underline once again that (|B3j) is the most general antisymmetric wave function. 

We express the functions $a and $5 as a superposition (we take V = 1 for simplicity) 



Vl-X 2 <Mri,r 2 ) = e ikiri+ik2r2 - e ^+^ , (B4) 



ki,k 2 

A$s(ri,r 2 ) = J2 ^.k, (e ikiri + ik2ra + e ikir2+ik2ri Y (B5) 

ki,k 2 ^ ' 

Therefore, the wave functions have the proper symmetry by construction. Setting either — or k 2 = 0> 

we obtain the singlet or triplet wave functions, respectively. Note that the binding energy for the triplet state vanishes 
for the s-wave pairing. We consider the Cooper problem with the wave function (|B3[) and Hamiltonian in the form 

As we already have said, the natural basis of the spin wave functions is spanned by the specific-spin wave functions. 
For this reason, we transform the wave function (|B3[) to the form 
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*(ri,r a ,(7i 1 CT 2 )= J2 [ (Ak 1)k2 e ikiri + ik2r2 +7 klik2 e ikir2 + ik2ri ) |1T)|21}- 

ki,k 2 

(7 kl ,k 2 e lkiri+ik2r2 + A klik2 e^ r -+*-^) |1 |)|2 ?)], (B6) 

where the new coefficients are given by 

A kljk2 = /?ki,k 2 + «ki,k 2 , (B7) 
7ki,k 2 = /3ki,k 2 — aki,k 2 - (B8) 

Following the standard but cumbersome procedure, we obtain the following two equations for those expansion 
coefficients 

1 E k (Ak + 7_ k+Q ^)^k 
V + 7 - k ' +Q ^ = N ^ + ^-2,-E ' (B9) 
, ,^ 1 S k (A k + 7- k+ Q^)V r (- k ^ +Q ^ f )(- k+Q A f) 

A k '+7- k ' +Q A f = R2Q2 ^ 2k , 2 — — . (BIO) 

2M + 2 Mred Z /- t 

For the above equations not to be contradictory, the matrix elements of the pairing potential have to be identical 

V ( _ k / + Q^p) (_k+Q^f») = Vuk. (Bll) 

If this condition is fulfilled, Eqs. (|B9|1 and (jBlOj) are equivalent. We introduce the potential in the standard form 
with a constant attraction 

v _ J -V , fork,k'e^*, m . 

Vkk '~\0, fork^WVk' <£W*. { > 

Now, the condition (|B11[) is fulfilled if and only if W* has the following property 

W * = -W* + Q — . (B13) 

Clearly, the region Wq is not a good choice for W* , but by making use of it we can construct a proper interaction 
region as follows 

/\ 777 

W* = W Q n(-W Q + Q— ). (B14) 
Such interaction region has a physical meaning, as we can rewrite it in the form 

vr= r k= kim 2 - k 2 mi i kl e Wl n WjJ A ka g ^ n Wl A q _ kj + k i ( B15) 
mi + m 2 

W t = {k|0<£ kffi <Kuj c }, i = l,2, (B16) 

meaning that the interaction takes place in the regions of k-space, for which not only both particles are at most 
hue Fermi level, but also if we exchanged their wavevectors (ki <->■ k 2 ), they are still at most Tiujc above their Fermi 
surface. 

The equation for the binding energy can be obtained from (|B9[) after summing up over k' 6 W* and dividing by 
^ kgM/ »(A k + 7_ k+ qAm). In effect, we have that 
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V E ¥W + ]fkL - 2u - E' (B17) 

This equation has the same form as the one obtained for the specific-spin state, the only difference being the 
appearance of the interaction region W* instead of Wq. From this very fact follows automatically the lower binding 
energy for the antisymmetric state (|B3I) . as W* C Wq. For the case of pair at rest (Q = 0), both states (specific-spin 
and antisymmetric) have equal binding energy, because then W* = Wq = 0. The numerical results show that the 
binding energy of the antisymmetric state decreases linearly with the increasing Q (see Fig. [5]). 



APPENDIX C: NARROW-BAND LIMIT OF THE ANDERSON LATTICE FOR U -» oo 



In this paper we consider the narrow-band limit of the Kondo-lattice state. Such limit has been considered earlier 25 ; 
here we provide a simple analytic argument. For the sake of simplicity we consider only h = case, but include also 
the hybrid pairing due to the Kondo-type coupling in the strong-correlation limit. The starting Hamiltonian with the 
term ~ V 2 /U included, has the form 



/ \ 2V 2 

E *m,ncL<w + E + v E + 4/*) - -^jj E (ci) 



m.n)(T 



e f + U 



where the first term represents the band energy of conduction electrons, the second the atomic energy of / electrons 
(N icr — fly fur), the third is the hybridization part, and the last the local pairing with ftt = (l/v^X/^cj^ — fl{ c \\)- 

The tilted operators exclude the double occupancies of / states, i.e. ji a = /i CT (l — flgfio) an d f\ a = — f\gfis)i U 
is the magnitude of the Hubbard / — / interaction. In the slave-boson saddle-point approximation this Hamiltonian 
representing the states in the lowest hybridized band reduces to 



'R = E^ata^ — s ■ = — al t al, a_k' i ak'Tj (C2) 

e/ + [/Z^ [(ek _ g/)2+ 4y 2 ]l/2[ (ek ,_ g/) 2 +4 y 2] i/2 k T ki M ^T. 

where 



£k- = i{e k + e> - [( £k - Iff + W 2 } 1/2 } (C3) 
and the if and V are the renormalized quantities e/ and V respectively 

Taking in the pairing part the k states close to the bare Fermi energy, we have that et — if ~ —if ~ V and under 
these circumstances the Hamiltonian (|C2[) reduces to the BCS form with Qk ff ps /k CT and a weakly k-dependent pairing 
potential, which is approximated by a constant Vq in main text (cf. Eq. (|29[) ). This approximation is justified also 
because the maximum pairing amplitude is achieved for ek = if, when the k-dependent ratio reduces to 4V 2 . The 
detailed estimate of E^ and of the pairing potential in terms of the effective Kondo temperature is more subtle^ 5 .. 
Roughly, the bandwidth of the heavy-quasiparticle band (of /-electrons) can be estimated as due to / — / hopping 
between the sites (ij) and is (V/ef) 2 tij — q(V/ef) 2 tij (since qV 2 = V 2 ). This is because it represents the three-step 
hopping sequence between neighboring /-states via conduction c-states: / — > c transition followed by a hopping in 
c-band and a subsequent c — > / deexcitation. One should underline that such simple form of the pairing part appears 
only if the main contribution to the hybridization V = Vim is of intraatomic [i = m) character. 
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